Accepted 2010 February 19 

Preprint typeset using MlgX style emulateapj v. 2/16/10 



O 



O 



> 

o 

(N 

P 

cn 
o 
o 



X 
J5 



DOPPLER SHIFT, INTENSITY, AND DENSITY OSCILLATIONS OBSERVED WITH THE EUV IMAGING 

SPECTROMETER ON HINODE 

John T. Mariska and K. Muglach^ 

space Science Division, Code 7673, Naval Research Laboratory, Washington, DC 20375; mariska@nrl.navy.mil 

Accepted 2010 February 19 

ABSTRACT 

Low- amplitude Doppler-shift oscillations have been observed in coronal emission lines in a number of active 
regions with the EUV Imaging Spectrometer (EIS) on the Hinode satellite. Both standing and propagating 
waves have been detected and many periods have been observed, but a clear picture of all the wave modes that 
might be associated with active regions has not yet emerged. In this study, we examine additional observations 
obtained with EIS in plage near an active region on 2007 August 22-23. We find Doppler-shift oscillations 
with amplitudes between 1 and 2 km s~^ in emission lines ranging from Fe XI 188.23 A, which is formed 
at logT = 6.07 to Fe XV 284.16 A, which is formed at logT = 6.32. Typical periods are near 10 minutes. 
We also observe intensity and density oscillations for some of the detected Doppler-shift oscillations. In the 
better-observed cases, the oscillations are consistent with upwardly propagating slow magnetoacoustic waves. 
Simultaneous observations of the Ca II H line with the Hinode Solar Optical Telescope Broadband Filter Imager 
show some evidence for 10-minute oscillations as well. 
Subject headings: Sun: corona — Sun: oscillations — Sun: UV radiation 



1. INTRODUCTION 

Because the outer regions of the solar atmosphere are 
threaded by a magnetic field, they can support a wide range 
of oscillatory phenomena. Theoretical aspects of these waves 
and oscillatio ns have been the subject of extensive investi- 
gations (e.g.. [Roberts et al.lll983l 11984 . These initial stud- 
ies were mainly driven by radio observations of short period 
oscillations (e.g.. [Rosen berg 1970). More recent detections 
of coronal oscillatory phenomena have resulted in consid- 
erable a dditional th eoretical work. Recent review s include 
iRobertsI (12000) and [Roberts & Nakariakovl (120031) . Obser- 
vational detections of coronal oscillatory phenomena include 
detections of spatial oscillations of coronal structures, which 
have been interpreted as fast kink mode MHD disturbances 
(e.g., Aschwanden et al. 1999); intensity oscillations, which 
have been i nterpreted a s propag ating slow magnetoacoustic 
waves (^e.g.. iDeForest & GurmanilT998l: iBerghmans & Clettd 
|1999|) : and Doppler shift oscillations, w hich have been in - 
terpreted as slow mode MHD waves (^e.g.. lWangetani2QQ2h . 
The interaction between the theoretical work and the grow- 
ing body of oscillation observations provides fertile ground 
for testing our understanding of the structure and dynamics of 
the corona. 

The EUV Imaging Spectrometer (EIS) on the Hinode satel- 
lite is an ex cellent tool fo r studying oscillatory phenomena in 
the corona. ICulhane eTaLl (|2007>) provides a detailed descrip- 
tion of EIS, and th e overall Hinode mission is described in 
iKosugi et al.l (|2007). Briefly, EIS produces stigmatic spectra 
in two 40 A wavelength bands centered at 195 and 270 A. Two 
slits (1'' and 2") provide line profiles, and two slots (40'' and 
266'0 provide monochromatic images. Moving a fine mirror 
mechanism allows EIS to build up spectroheliograms in se- 
lected emission lines by rastering a region of interest. With 
typical exposure times of 30 to 90 s, however, it can take con- 
siderable time to construct an image. Higher time cadences 
can be achieved by keeping the EIS slit or slot fixed on the 
Sun and making repeated exposures. This sit-and-stare mode 

^ Also at Artep Inc., EUicott City, MD 21042 



is ideal for searching for oscillatory phenomena. 

EIS Doppler shift data have already been used for 
a number of investi g ations of oscillatory phenomena. 
IVan Do orsselaere et al.l (12008 ^ have detected kin k mode 
MHD oscillations with a period near 5 minutes. Mari ska et al.l 
(2008) observed damped slow magnetoacoustic standing 
waves with periods of about 35 minutes. IWang et aH (2009a) 
have detected slow mode magnetoacoustic waves with 5 
minute periods propagating upwar d from the c hromos phere 
to the corona in an active region. IWang et aH (l2009b|) have 
also observed propagating slow magnetoacoustic waves with 
periods of 12 to 25 minutes in a large fan structure associated 
with an active region. Analysis of oscillatory data with EIS 
is still just beginning. The amplitudes observed have all been 
very small — typically 1 to 2 km s~^. Thus a clear picture of 
the nature of the low-amplitude coronal oscillations has yet 
to emerge. In this paper, we add to that picture by analyzing 
portions of an EIS sit-and-stare active region observation that 
shows evidence for Doppler-shift oscillations in a number of 
EUV emission lines. We also use data from the Hinode Solar 
Optical Telescope (SOT) to relate the phenomena observed in 
the corona with EIS to chromospheric features and the mag- 
netic field. 

2. OBSERVATIONS 

The observations discussed in this paper were taken in an 
area of enhanced EUV and soft X-ray emission just west of 
NOAA active region 10969. The complete EIS data set con- 
sists of a set of 256'' x 256" context spectroheliograms in 20 
spectral windows, followed by 7.4 h of sit-and-stare observa- 
tions in 20 spectral windows, and finally a second set of con- 
text spectroheliograms. All the EIS sit-and-stare data were 
processed using software provided by the EIS team. This pro- 
cessing removes detector bias and dark current, hot pixels, 
dusty pixels, and cosmic rays, and then applies a calibration 
based on the prelaunch absolute calibration. The result is a 
set of intensities in ergs cm~^ s~^ sr~^ A~^ The data were 
also corrected for the EIS slit tilt and the orbital variation in 
the line centroids. For emission lines in selected wavelength 
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Figure 1. Context EIS spectroheliograms in 9 wavelength windows obtained from 13:33:46 to 17:56:57 UT on 2007 August 22 (left) and from 01:55:43 to 
06:18:53 UT on 2007 August 23 (right). The location of the EIS slit for the sit-and-stare observation is marked with the vertical line. The horizontal lines along 
the slit mark the location of the oscillation data analyzed in this paper. The post-flare loops discussed in the text are the bright features east and south of the 
marked region in the three bottom right panels. 



windows, the sit-and-stare data were fitted with Gaussian fine 
profiles plus a background, providing the total intensity, loca- 
tion of the line center, and the width. 

Figure [T] shows the pre- and post- sit- and- stare spectroheli- 
ograms and provides a more detailed view of the area on the 
Sun covered by the sit-and-stare observations. The spectrohe- 
liograms were obtained with the V EIS slit using an exposure 
time of 60 s. Note that the two sets of spectroheliograms have 
small differences in the pointing. Examination of the two sets 
of images shows that considerable evolution has taken place 
in the observed region between the times each set was taken. 
In particular, new loop structures have developed at locations 
to the south of the bright core of the emitting region, suggest- 
ing that some flare-like heating may have taken place between 
the two sets of context observations. Space Weather Predic- 
tion Center data show that there was one B1.2 class flare be- 
tween 07:50 and 08: 10 UT on 2007 August 22 at an unknown 
location. Since other flares occurred in this active region, it 
is likely that this one is also associated with it. There were, 
however, no events recorded during the time of the EIS obser- 
vations. 

Both magnetogram data from the SOHO Michelson 
Doppler Imager (MDI) and Hinode Solar Optical Telescope 
(SOT) and coronal images taken in the 171 A (Fe IX/X) fil- 
ter with the Transition Region and Coronal Explorer (TRACE) were 
also examined. The magnetograms, one of which is shown in Fig- 
ure|2]with the location of the EIS slit for the sit-and-stare observation 
indicated, show an extended bipolar plage region without any visible 
sunspots. As can be seen in the magnetogram, the EIS slit crosses 
two regions of positive (white) polarity. The weaker one, which has 
a size of about 6''and is marked with an arrow, is the focus of this 
analysis. 

The TRACE images, one of which is shown in Figure[3] show that 
the western part of the region near the right edge of the box that out- 
lines the edges of the second EIS spectroheliogram exhibits fan-like 
coronal structures extending radially out in all directions, while the 
eastern part consists of shorter loop structures that connect the oppo- 
site polarities. The full set of TRACE images that cover most of the 
time period covered by the context and sit-and-stare observations are 
included as a movie. The movie shows generally quiescent behav- 



ior in the fan-like structures, but multiple brightenings in the eastern 
loop system. The largest of these took place at around 01:31 UT on 
August 23 and shows a filament eruption, which resulted in a system 
of post-flare loops that are visible in the Fe XV 284.16 A and Fe XVI 
262.98 A spectroheliograms in the right panels of Figure[T] Some of 
the loops in the eastern part cross the EIS slit location for the sit-and- 
stare observations, but they are for the most part north of the area we 
focus on. 

We have also examined images taken with the Extreme Ultravi- 
olet Imaging Telescope (EUVI, Wuelser et al. 2004) on the Ahead 
(A) and Behind (B) spacecraft of the Solar Terrestrial Relations Ob- 
servatory (STEREO). EUVI observes the entire solar disk and the 
corona up to 1.4 R© with a pixel size of 1.59'^ On 2007 August 22, 
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Figure 2. A portion of an MDI magnetogram taken on 2007 August 22 at 
14:27 UT showing the area covered by the EIS spectroheliograms and the 
approximate location of the EIS slit during the sit-and-stare observation. The 
horizontal lines along the slit mark the location of the oscillation data ana- 
lyzed in this paper, and the area of weak magnetic flux that is the focus of 
this analysis is marked with an arrow. The image has been scaled so that the 
range of magnetic fluxes displayed is ±100 Gauss. 
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Figure 3. TRACE 171 A (Fe IX/X) image taken on 2007 August 22 
at 20:01:19 UT. The box shows the location of the EIS spectroheli- 
ograms taken before and after the sit- and- stare observations, and the 
vertical line with the box shows the location of the sit-and stare ob- 
servation. The horizontal lines along the slit mark the location of the 
oscillation data analyzed in this paper. 



the separation angle of STEREO with Earth was 14.985° for space- 
craft A and 11.442° for spacecraft B. Our target region was rather 
near the limb as seen from spacecraft A. EUVI/B, however, provided 
continuous on-disk images between 2007 August 22, 11:00 UT and 
2007 August 23, 11:30 UT, which allowed us to study the develop- 
ment of the coronal structures of the region during the time of the 
EIS observations. 

We produced movies in all four EUVI channels (171 A, 195 A, 
284 A, and 304 A), with 171 A having the highest cadence of 2.5 
min. The EUVI movies generally show the same behavior as that 
shown in the TRACE data. EUVI captured two larger events during 
this period, the first one started at around 12:24 UT on August 22, 
showing multiple brightenings of the eastern loop system although 
it is not clear if a CME had been launched. The second one took 
place at around 01:31 UT on August 23, involving the western part of 
the region and shows a filament eruption with a system of post-flare 
loops evolving at the time the second EIS context spectroheliogram 
was taken. 

SOT obtained Call H observations at a one-minute cadence during 
the time interval of the sit-and-stare observations that are the focus of 
this study. To coalign those observations we use an EIS O VI 184.12 
A spectroheliogram, MDI magnetograms, SOT Spectro-Polarimeter 
(SP) magnetograms, and SOT Ca II H line filtergrams. The EIS 
O VI 184.12 A spectroheliogram was part of the scan shown in Fig- 
ure[T](left), taken several hours before the sit-and-stare observations. 
Coalignment was carried out by rescaling the images, using cross- 
correlation techniques to overlap the selected subfields and finally 
blinking the images to ensure minimal offsets. 

Figure |4] displays the O VI 184.12 A spectroheliogram in the top 
left corner, with the dark vertical line representing the slit position of 
the sit-and-stare observations. The top right image is from a full-disk 
MDI magnetogram taken at 14:27 UT. The three small structures in 
the lower left corner of the magnetogram were used to coalign the 
O VI image, which shows bright emission structures at the same lo- 
cations. In the next step we use SP data to coalign with the MDI mag- 
netogram. The SP was scanning the target region on 2007 August 22 
between 18:03 UT and 19:01 UT, at the same time the Ca II H image 
sequence was taken. The lower right image in Figure |4] gives the ap- 
parent longitudinal field strength as derived from the Stokes spectra. 
Comparing the two magnetograms we can see that the overall struc- 



ture of the bi-polar region has hardly evolved over time and that the 
SP data show considerable additional fine structure in the magnetic 
field. We also note that in the time between the MDI magnetogram 
and SP magnetogram flux cancellation has taken place at the polar- 
ity inversion line of the bi-polar region (at around x = —610^'). The 
eruption observed a few hours later in EUVI and TRACE was proba- 
bly triggered by this flux cancellation. 

The target area for the oscillation study is in the black circle, a 
small magnetic region of positive polarity. In the final coalignment 
step the SOT Ca II H images were coaligned with the SP scan. The 
bottom left image in Figure|4]shows the temporal average of the time 
sequence taken between 19:50 UT and 21:40 UT. The black circle 
encompasses a small network bright point that overlaps with the EIS 
slit and the region averaged along the slit in the analysis that is the 
focus of this paper. As we have already noted, the EIS spectrohe- 
liograms taken before and after the sit-and-stare observation show 
that considerable structural evolution has taken place. In the Fe XV 
264.78 A spectroheliogram taken after the sit-and-stare observation, 
there is significant emission within the 6'' region of the EIS slit that 
is the focus of this paper, suggesting that the location is near one end 
of one of the coronal loops seen in the lower three panels in the right 
side of Figure [T] On the other hand the images in the EIS spectro- 
heliograms taken immediately before the sit-and-stare observation 
show much less evidence for a loop rooted at that location. Thus we 
can only conclude that the bright feature we observe in the Ca II H 
line may be the footpoint of a coronal loop. We do note, however, 
that there are no other strong magnetic concentrations close to the 
location of the region that is the focus of this study. The ones to 
the south-east are about 10'' away — more than our co-alignment er- 
rors. We believe that loop footpoints will be anchored in strong flux 
concentrations. 

The basic building block for the sit-and-stare observation was EIS 
Study ID 56. This study takes 50 30 s exposures in 20 spectral 
windows with a solar j-position coverage of 400''. The entire sit- 
and-stare observation consisted of four executions of the study, with 
each execution invoked with a repetition count of four. Thus, the full 
sit-and-stare data set consists of 16 repetitions of the basic building 
block. There was a brief delay between each invocation of the study, 
leading to gaps in the resulting time series. 

Over both orbital periods and shorter time intervals, the Hinode 
pointing fluctuates in both the solar x- and j-directions by up to about 
3" peak-to-peak. This is due to both spacecraft pointing variations 
and changes in the location of EIS relative to the location of the 
Sun sensors on the spacecraft. Studies of the latter variations have 
shown that they are well-correlated with similar variations observed 
in data from the Hinode X-Ray Telescope (XRT). We have therefore 
used a modified version of the software provided by the XRT team to 
compute the average j-position of each pixel along the slit over each 
exposure and then interpolated the fitted centroid positions onto a 
uniform j-position grid. 

It is not, of course, possible to correct the sit-and-stare observa- 
tions for fluctuations in the x-position on the Sun. Plots of the x- 
position pointing fluctuations show that they tend to be smoother 
than the j-direction fluctuations, and that they are dominated by the 
orbital period. If periodic Doppler shifts are present only in very 
small structures, we would thus expect the signal to show a mod- 
ulation with the orbital period. Larger structures, > 3" in the x- 
direction, that display coherent Doppler shifts should not be affected 
by the spacecraft x-position pointing variations. 

Over much of the EIS slit during the sit-and-stare observation 
there is no evidence for interesting dynamical behavior in the 
Doppler shift observations. The portion of the slit that covers the 
brighter core area of the region does, however, show evidence for 
periodic changes in the Doppler shifts. Figure [5] shows the measured 
Doppler-shift data in four emission lines over this j-position range 
as a function of time. The gaps between each invocation of the study 
appear as wider pixels near 20, 22, and UT. Note that the data are 
also affected by passage of the spacecraft through the South Atlantic 
Anomaly (SAA). A small region of compromised data appears near 
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Figure 4. Coalignment of EIS and SOT. Top left: EIS spectroheliogram taken in O VI 184.12 A at the same time as the spectroheHograms shown in Figure [T] 
the dark vertical line near the center marks the slit location for the sit-and-stare observations. Top right: MDI full disk magnetogram taken at 14:27 UT. Bottom 
right: Hinode SP magnetogram, taken between 18:03 UT and 19:01 UT, at the same time the Ca 11 H time sequence was taken. Bottom left: temporal average of 
the time sequence in Ca 11 H taken between 19:50 and 21:40 UT. The black circles in the bottom images encompass a small network bright point that overlaps 
with the EIS slit. 



19:30 UT, and major SAA passages are evident near 21:00 and 22:45 
UT. 

The display shows clear evidence for periodic fluctuations in all 
the emission lines at j-positions of roughly —170'' to —180'' during 
the first hour shown in the plots. Periodic fluctuations are also visi- 
ble over j-locations centered near —210". Note that there is some 
evidence of longer period fluctuations, for example in the Fe XV 
284.16 A emission line. These fluctuations are probably the result 
of the correction for orbital line centroid shifts not fully removing 
those variations. In the remainder of this paper, we focus on the area 
near —210" that shows evidence for oscillatory phenomena. 

3. ANALYSIS 

3.1. Doppler Shift Oscillations 

As we noted earlier, solar j-positions between —200" and —215" 
show considerable oscillatory behavior, particularly in the set of 
data taken beginning at 20:07:01 UT. Figure [6] shows the averaged 
Doppler shift data over the time period of this set of observations 



for, from top to bottom, Fe XI 188.23 A, Fe XII 195.12 A, Fe XIII 
202.04 A, Fe XIV 274.20 A, and Fe XV 284. 16 A. The Doppler shifts 
have been averaged over the 16 detector rows from —200" to —215". 
Data from 172 to 180 minutes were taken during SAA passage and 
have been removed from the plot. 

The Doppler shift data over this portion of the EIS slit show clear 
evidence for low-amplitude, roughly 2-4 km s~^, oscillatory behav- 
ior with a period near 10 minutes. For some of the time period, 
particularly after 180 minutes, there appears to be a clear trend for 
the oscillations to display increasing amplitude as a function of in- 
creasing temperature of line formation. 

Because there is a significant gap in the Doppler shift data, nei- 
ther Fourier time series analysis nor wavelet analysis is appro- 
priate. Instead, we examine the time series by calculating peri- 
odograms using the approach outlined in Home & Baliunas ( 198^ 
and Press & Rybicki ( 1989). Figure [T] shows the periodograms cal- 
culated from the Doppler shift data shown in Figure [6] Also plotted 
on the figure are the 99% and 95% significance levels. All but one 
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Figure 5. EIS sit-and-stare Doppler-shift data in four emission lines covering 
the entire sit-and-stare observation. The Doppler shift in each window has 
been adjusted so that the zero value is the average over the window. The 
maximum and minimum values plotted in each case are +20 and —20 km 
s~\ respectively. This study focuses on the spatial range from —200" to 
-215" in the time interval from 20:07 to 21:53 UT. 
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Figure 6. Doppler shift data averaged over the 16 detector rows from —200" 
to —215". The emission lines shown from top to bottom are Fe XI 188.23 A, 
Fe Xll 195.12 A, Fe Xlll 202.04 A, Fe XIV 274.20 A, and Fe XV 284.16 A. 
Time is measured in minutes from the start of the sit-and-stare data at 
18:13:06 UT. The zero value for Doppler shifts is set to the average Doppler 
shift in each line and is shown using the horizontal lines. Data for each emis- 
sion line are displaced from the next set by 4 km s~^ 



of the time series show a peak in the periodogram at the 99% con- 
fidence level, and the largest peak in the Fe XV 284.16 A emission 
line data is at the 95% confidence level. Table [T] lists the period of 
the most significant peak in each of the panels shown in the figure. 

To estimate the amplitude of the oscillations, we detrend the 
Doppler shift time series by subtracting a background computed us- 
ing averaged data in a 10-minute window centered on each data point 
and then compute the rms amplitude as the standard deviation of the 
mean. For a sine wave, the peak velocity is the rms value multiplied 



Table 1 

Periods and Amplitudes Detected in Doppler Shift and Intensity Data 





Wavelength 


Logr 


Pd 


5v 


Pi 


5I/I 


Ion 


(A) 


(K) 


(minutes) 


(km s-i) 


(minutes) 


(%) 


Fe XI 


188.23 


6.07 


10.0 


1.2 


11.8 


1.1 


Fe XII 


195.12 


6.11 


10.1 


1.1 


11.4 


0.9 


Fe XIII 


202.04 


6.20 


9.1 


1.3 


11.2 


1.4 


Fe XIV 


274.20 


6.28 


9.1 


1.3 


11.4 


2.0 
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284.16 


6.32 


9.0 


1.4 
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2.4 
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5 10 15 

Period (min) 

Figure 7. Periodograms computed for the Doppler shift data shown in Fig- 
ure [51 The solid horizontal line on each plot indicates the power equivalent 
to a calculated false alarm probability of 1% and the dashed line is for a false 
alarm probability of 5%. 

by V2. These peak velocity values are also listed in the table in the 
dv column. Visual inspection of the data suggests that the numbers 
in the table are smaller than what might be obtained by fitting the 
data. The numbers do confirm the impression that the oscillation 
amplitude increases with increasing temperature of line formation. 

It is clear from the Doppler shift plots in Figure [6l that the oscilla- 
tions are not always present. Instead, they appear for a few periods 
and then disappear. To understand better this behavior, we have fitted 
the time intervals where oscillations are obvious with a combination 
of a damped sine wave and a polynomial background. Thus, for 
each time period where oscillations are present we assume that the 
data can be fitted with a function of the form 



v(^) = Ao sin(co^ + 0) exp(-AO ^B{t) 



(1) 



where 

B{t) = bQ + bit + b2t^ + Z?3^^ + ■ • • (2) 

is the trend in the background data. Time is measured from an ini- 
tial time tQ, which is different for each set of oscillations we fit. The 
fits were carried out using Levenberg-Marquardt least- squares min- 
imization (Bevington 1969). Generally only two terms in the back- 
ground polynomial were necessary. 

Figure [S] shows the results of this fitting for the Doppler shift data 
beginning 142.5 minutes after the start of the sit-and-stare observa- 
tion. All the fits show roughly the same amplitudes, periods, and 
phases. Emission lines formed at the higher temperatures (Fe XIV 
and Fe XV) show clear evidence for more than one full oscillation 
period. At lower temperatures, the oscillatory signal damps much 
more rapidly. 

Table [3 lists the amplitudes, Aq, periods, P, phases, 0, and the 
inverse of the decay rate. A, that result from fitting all the time pe- 
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Figure 8. Decaying sine wave fits to the Doppler shift data beginning roughly 
142 minutes after the start of the sit-and-stare observation. The plotted data 
has had the polynomial background removed. Vertical dashed lines show the 
time range used in the fitting. 



Table 2 

Doppler Shift Oscillation Properties 



to 




Ao 


P 






(min) 


Ion 


(km s-i) 


(min) 


(rad) 


(min) 


113.9 


Fe XI 


0.9 


8.8 


1.8 


19.0 


113.9 


Fe Xll 


0.7 


9.2 


2.3 


-216 


113.9 


Fe Xlll 


0.7 


8.8 


2.3 


-73.3 


113.9 


Fe XIV 


0.3 


10.7 


3.9 


-10.9 


113.9 


Fe XV 


0.5 


9.5 


2.9 


-24.5 


142.5 


Fe XI 


2.4 


10.4 


1.0 


5.1 


142.5 


Fe XII 


2.0 


8.4 


0.9 


6.6 


142.5 


Fe XIII 


1.6 


8.2 


1.2 


15.0 


142.5 


Fe XIV 


2.0 


8.7 


1.9 


23.9 


142.5 


Fe XV 


2.2 


8.9 


2.1 


23.3 


190.0 


Fe XI 


0.3 


8.8 


4.7 


-411 


190.0 


Fe XII 


0.8 


7.5 


3.5 


36.8 


190.0 


Fe XIII 


1.2 


7.3 


3.1 


36.3 


190.0 


Fe XIV 


1.5 


7.3 


3.2 


43.6 


190.0 


Fe XV 
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riods in the data for which a reasonable fit to the Doppler shift data 
could be obtained. The periods are consistent with the results of the 
periodogram analysis for the entire time interval. Generally, the am- 
plitudes are larger than the 3v values show in Table[T] Note that some 
of the fits show negative decay times, indicating that some of the os- 
cillations show a tendency to grow with time. In these cases, this 
is not followed by a decay, but rather a rapid loss of the oscillatory 
signal. 

All the fits use the start times listed in the table and thus the 
phase values for each time interval can not be directly compared. 
When the phases are adjusted to a common start time, the values do 
not agree. Thus while the periods are similar for each time interval 
in which oscillations are observed, it appears that each event is being 
independently excited. 

Examining the amplitude of the oscillations in each time interval 



as a function of the temperature of line formation shows no clear 
trend. The data set starting at 113.9 minutes seems to show evidence 
for a decrease in amplitude with increasing temperature, while the 
data set beginning at 190.0 minutes shows the opposite trend. Simi- 
larly, there is no clear trend in the periods of the oscillations in each 
data set as a function of temperature. 

3.2. Intensity Oscillations 

If the observed Doppler shift oscillations are acoustic in nature, 
then they should also be visible in the intensity data. For a linear 
sound wave v = csSp/p, where v is the amplitude of the wave, cs 
is the sound speed, and Sp is the density perturbation on the back- 
ground density p. Taking an amplitude of 2 km s~^, yields values of 
5p/p of around 1%. Since the intensity fluctuation, 51 /I, is propor- 
tional to 25p/p, we expect only about a 2% fluctuation in the mea- 
sured intensity. This number could of course increase if the actual 
velocity is much larger due to a large difference between the line- 
of- sight and the direction of the coronal structure being measured. 
Figure |9] shows the measured intensity data averaged over the same 
locations as the Doppler shift data shown in Figure [6] The data show 
little or no evidence for oscillations with the periods measured in the 
Doppler shift data. This is borne out by a periodogram analysis of 
the time series in the figure, which show no significant peaks. 

If, however, we detrend the data by subtracting the gradually 
evolving background signal, there is some evidence for an oscilla- 
tory signal. Figure [To] shows the data in Figure|9]with a background 
consisting of a 10-minute average of the data centered on each data 
point subtracted. All the emission lines show some evidence for an 
oscillatory signal, with the Fe XIII 202.04 A emission line being the 
most obvious. 

Figure [TT] shows periodograms constructed for the emission lines 
shown in Figure [Tol Each periodogram shows a significant peak. 
The periods for the strongest peak in each periodogram are listed in 
Table [U The periods are generally consistent with those determined 
from the Doppler shift data. Also listed in the table is an estimate 
of the intensity fluctuation in each emission line. This was obtained 
by computing the standard deviation of the detrended intensity in 
the time series for each emission line and then dividing the result by 
the average intensity. The values are roughly consistent with those 
expected based on the 5v estimates listed in Table [T] 

While the oscillatory signal is much less strong in the detrended 
intensity data than in the Doppler shift data, it is possible to fit some 
of the data in roughly the same time intervals that were used for the 
results listed in Tabled Table [3] shows the resulting fit information. 
Note that intensity oscillation fits were not possible for all the lines 
for which the Doppler shift data could be fitted. Figure [12] shows an 




120 140 160 180 200 220 
Time (min) 

Figure 9. Normalized intensity data averaged over the 16 detector rows from 
—200" to —215". The emission lines shown from top to bottom are Fe XI 
188.23 A, Fe XII 195.12 A, Fe XIII 202.04 A, Fe XIV 274.20 A, and Fe XV 
284.16 A. Time is measured in minutes from the start of the sit-and-stare data 
at 18:13:06 UT. Data for each emission line are displaced from the next set 
by 0.2. 
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Figure 10. Detrended intensity data averaged over the 16 detector rows from 
—200" to —215". The emission lines shown from top to bottom are Fe XI 
188.23 A, Fe XII 195.12 A, Fe XIII 202.04 A, Fe XIV 274.20 A, and Fe XV 
284.16 A. Time is measured in minutes from the start of the sit-and-stare data 
at 18:13:06 UT. Data for each emission line are displaced from the next set 
by 100 ergs cm~^ s~^ sr~^. 
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Figure 11. Periodograms computed for the detrended intensity data shown in 
normalized form in Figure |9] The solid horizontal line on each plot indicates 
the power equivalent to a calculated false alarm probability of 1% and the 
dashed line is for a false alarm probability of 5%. 

example of the fits to the detrended intensity data. For these plots 
the intensity has been converted to a residual intensity expressed 
in % by taking the difference between each data point and the 10- 
minute average and dividing it by the 10-minute average. The plot- 
ted points have also had the polynomial background subtracted.To 
facilitate comparisons with the Doppler shift data we have also plot- 
ted as green curves on the panels for the Fe XIV and Fe XV intensity 
data the fitting results for the Fe XIV and Fe XV Doppler shift data. 

Comparison of the fits in Figure [T2l with those in Figure [8] shows 
some similarities and many differences between the two data sets. 
For the lines with a reasonably strong signal, the periods measured 
for the residual intensity data are similar to those measured for 
the Doppler shifts. In addition, the residual intensities show the 
same trend to larger amplitudes as the temperature of line forma- 
tion increases. The intensity oscillations clearly start earlier than the 
Doppler shift oscillations. Moreover, while the Doppler shift oscilla- 
tions are damped, the intensity oscillations appear to grow with time 
over the same interval. This appears to be the case for the other time 
intervals as well. Fitting the detrended intensity signal is more chal- 
lenging than fitting the Doppler shift signal. Also, implicit in the 




Figure 12. Decaying sine wave fits to the detrended intensity data beginning 
roughly 138 minutes after the start of the sit-and-stare observation. The in- 
tensities have been converted to residual intensities by taking the difference 
between the intensity at each data point and the running mean intensity and 
dividing it by the running mean. The plotted data has also had the polyno- 
mial background used in the fit removed. Vertical dashed lines show the time 
range used in the fitting. Also plotted as green curves on the panels for the 
Fe XIV and Fe XV intensity data are the fitting results for the Fe XIV and 
Fe XV Doppler shift data. Those plots show —v{t). 
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fitting is the idea that a damped sine wave can fully represent what is 
probably a more complex signal. Thus, we are reluctant to read too 
much into this growth in the intensity until we can determine from 
additional data sets if it is a common phenomenon. 

An important factor in determining the nature of the oscillations 
is the phase difference between the Doppler shift signal and the in- 
tensity signal. Comparing the phases of the fits listed in Table|2]with 
those in Table[3]is difficult because the periods are not identical, but, 
since the periods are close, the differences do not significantly alter 
any conclusions that we might draw. To facilitate this comparison we 
have plotted as green curves on the panels for the Fe XIV and Fe XY 
intensity data in Figure[T2]the fitting results for the Fe XIV and Fe XV 
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Doppler shift data. For the Doppler shift data, we plot — v(^). The 
curves show that for these two ions, the intensity variations are close 
to 180° out of phase with the Doppler shift variations. Since we de- 
fine the Doppler shift as c3X/X, this means that the peak intensity 
corresponds to a blueshift, indicating an upward propagating wave. 
For the other two time intervals, the situation is more ambiguous. 
Examination of the tables shows that in many cases, the periods are 
significantly different for the same Doppler shift and intensity data in 
the same line. In those cases where the periods are close (e.g., Fe XII 
at ^0 = 1 13.9 minutes and Fe XV at ?o = 190.0 minutes), examination 
of the plots similar to Figures [Sl and [T2l shows the same 180° phase 
shift, again indicating upwardly propagating oscillations. 

Even for the cases where the periods are close, the agreements in 
the phases are only approximate. For the Fe XIV and Fe XV intensity 
and Doppler shift fits shown in Figure [TJl the intensity oscillation 
leads the Doppler shift oscillation by a small fraction of a period. 
For both the Fe XII data at = 113.9 minutes and the Fe XV data 
at t{) = 190.0 minutes, the Doppler shift oscillation leads the inten- 
sity oscillation by a fraction of a period. In both cases, this differ- 
ence is les s than the 1/4 p e riod expected for a standing-mode MHD 
wave (Sak urai et al .1120021) . IWang et aP (^009a) observed propagat- 
ing waves with periods in the four to six minute range in EIS active 
region observations. They noted that for most cases the Doppler shift 
and intensity oscillations were nearly in phase. In the cases where 
there was a difference, the phase of the intensity was earlier than the 
Doppler shift, as is the case for the data shown in Figures [8l and [T2l 
Theoretical modeling of propagating slow waves with periods near 
five minutes show that thermal conducti on can produce p hase shifts 
between the intensity and Doppler shifts (lOwen et al.l2009 ). Further 
study of EIS data sets where both the Doppler shift and intensity can 
be fitted could provide valuable constraints on these models. 

3.3. Density Oscillations 

The electron density is one factor in determining the Alfven speed 
in the oscillating plasma. Moreover, for magnetoacoustic fluctua- 
tions, we expect the density to also oscillate. Thus a direct measure- 
ment can aid in disentangling the nature of the oscillations. The 
sit-and-stare observations included density- sensitive line pairs of 
Fe XII (186.88 A/195.12 A) and Fe XIII (203.83 A/202.04 A). Using 
data from version 5.2 of the CHIANTI database (Landi et al. 2006; 
iDere et alT l997). we computed the electron density at each time for 
the row-averaged data. Fo r Fe XII, C HIANTI uses energy levels and 
radiative decay rates from Del Zanna & Mason (2005), electron col- 
lision strengths from Storey et al. (2005), and proton collision rate 
coefficients from Landman (1978). Fo r Fe XIII, CHIANT I uses en- 
ergy levels from Penn & Kuhn ( 19941). lJupen et al.l (Il993h . a nd ver- 
sion 1 .0 of the NIST database; radiative decay rates from 'Young* 
(120041) . electron collision strengths from Gupta & Taval ( 1998), and 
proton collision rate coefficients from Landman ( 1975). These diag- 
nostics are discussed in detail in Young et al. (2009). 

Figure [13] shows the derived electron densities as a function of 
time for the same time interval shown in Figures [8] and [12] Both 
sets of derived densities show the same overall time behavior, but 
the absolute values differ by nearly a factor of three. Differences 
between these diagnostics have been noted before (e.g.. Young et al. 
^009), and are thought to be due to issues with the atomic data. It 
is not yet clear which of the values should be considered the most 
reliable. 

Neither of the density time series shown in the figure display any 
evidence for the oscillations detected in the detrended intensity data. 
If we smooth the density time series with a 10-minute running mean, 
there is some evidence for oscillatory behavior over the time range 
beginning at 190 minutes. Decaying sine wave fits to that region are 
show in Figure [14] For the Fe XII time series the fit has an amplitude 
of 1.9 X 10^ cm~^, a period of 13.5 minutes, a decay time of —46.4 
minutes, and a phase of 3.5 radians, all consistent with the values 
listed for the Fe XII data in Table [3] for this time interval. For the 
Fe XIII time series the fit has an amplitude of 5.1 x 10^ cm~^, a 
period of 10.9 minutes, a decay time of —25.1 minutes, and a phase 
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Figure 13. Electron density determined using the Fe XII 186.88 A/195.12 A 
ratio (top) and the Fe XIII 203.83 A/202.04 A ratio (bottom) for the data 
beginning roughly 138 minutes after the start of the sit-and-stare observation. 
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Figure 14. Decaying sine wave fits to the smoothed electron density de- 
termined using the Fe XII 186.88 A/195.12 A ratio (top) and the Fe XIII 
203.83 A/202.04 A ratio (bottom) for the data beginning 190 minutes after 
the start of the sit-and-stare observation along with. The plotted data has had 
the polynomial background removed. Vertical dashed lines show the time 
range used in the fitting. 



of 1.9 radians. With the exception of the phase, these values are 
generally consistent with the Fe XIII data in Table [3] for this time 
interval. The amplitudes for the Fe XII and Fe XIII fits are 0.5% and 
0.4% of the average density in the time interval. Since the observed 
intensity fluctuations are only about 1%, though, and we expect 51 /I 
to be proportional to 25p/p, they are consistent with the observed 
intensity oscillations. 

3.4. Underlying Chromo spheric Behavior 

As we pointed out earlier, SOT obtained a time sequence of Ca II 
H images that is co-temporal with the EIS sit-and-stare observations 
starting at 19:50 UT and ending at 21 :40 UT, with a constant cadence 
of 60 s. To further use this data, we applied the standard reduction 
procedure provided by the SOT team that is available in SolarSoft. 
The images of the Ca II H sequence were then carefully aligned using 
Fourier cross-correlation techniques to remove residual jitter and the 
drift of the SOT correlation tracker. 

As can be seen in the lower left panel of Figure [4] the EIS slit cov- 
ers the network bright point, which is the focus of the chromospheric 
analysis. Considering the accuracy of the coalignment and the spa- 
tial averaging applied to the EIS data, we also average the Ca II H 
signal. 

Figure [15] shows the time history of the Ca II H-line intensity for 
three different sized spatial areas centered on the feature shown in 
Figure [4] The sizes of the regions are given in SOT pixels, which 
are 0.10896^' in size. As expected for chromospheric lines, all three 
averaged data sets show evidence for intensity oscillations with a 



SOLAR DOPPLER SHIFT, INTENSITY, AND DENSITY OSCILLATIONS 



9 



100 



c 

CD 




140 160 
Time (min) 

Figure 15. Ca II H intensity data averaged over different spatial areas cen- 
tered on the magnetic feature highlighted in Figure |4] Each data set has 
had the mean value subtracted. Spatial areas are given in SOT pixels, 
which are O.I0896"in size, giving areas of 2.18" x 2.18", 3.77" x 3.77", and 
4.36" X 4.36", respectively. 



1 1 1 1 1 1 1 1 1 1 

- 40x40 









1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- 30x30 


: J 








5 10 
Period (min) 

Figure 16. Periodograms computed for the detrended Ca II H intensity data. 
The size of the area averaged over in pixels in indicated in each plot. The solid 
horizontal line on each plot indicates the power equivalent to a calculated 
false alarm probability of 1% and the dashed line is for a false alarm proba- 
bility of 5%. Spatial areas are given in SOT pixels, which are 0.10896"in size, 
giving areas of 2.18" x 2.18", 3.77" x 3.77", and 4.36" x 4.36", respectively. 



period near 5 minutes. The data, however, are quite noisy and pe- 
riodograms constructed from them show no significant peaks. De- 
trending the data, however, results in periodograms with significant 
peaks. Figure [161 shows periodograms for the three data sets with 
each detrended by subtracting a 9-minute running mean from each 
data point. The periodograms show clear evidence for oscillations 
near 5 minutes. 

There is also some evidence for power at periods between 9 and 
10 minutes, but with more than a 5% false alarm probability. If we 
instead detrend the data with an 11 -minute running mean, then the 
peak between 9 and 10 minutes becomes more prominent and is at or 
above the 5% false alarm probability level for the 30 pixel x 30 pixel 
and 20 pixel x 20 pixel data sets. Wavelet analysis of the data sets 
detrended with both a 9- and 11 -minute running mean shows signif- 
icant signal near 9 minutes. 

Note that the data in Figure [15] shows three larger peaks in all 



three data sets. The first peak, near 110 minutes, occurs before the 
beginning of the EIS data plots in Figures[6land[T0l The other two, at 
roughly 140 and 170 minutes, come just before significant oscillatory 
signal is observed in the EIS Doppler shift and detrended intensity 
data. It is tempting to suggest that these enhancements correspond to 
chromospheric events that resulted in the oscillations observed with 
EIS. 

In principle the EIS He II 256.32 A data can bridge the gap in 
temperature between the SOT Ca II data and the Fe lines formed at 
higher temperatures that are the main focus of this study. In practice 
the He II data are challenging to analyze. The line is closely blended 
with a Si X line at 256.37 A along with a smaller contr ibution from 
Fe X and Fe XIII ions at slightly longer wavelengths dBrown et al.1 
12008 1). In an effort to see if a connection can be made, we have 
made two-component Gaussian fits to the row-averaged He II data. 
Periodograms of the resulting Doppler shift data show no significant 
periods. Periodograms of the He II fitted intensities detrended with a 
10-minute running mean show no peaks at the 1% false alarm prob- 
ably level and one peak with a period between 12 and 13 minutes at 
the 5% false alarm probability level. Examining the detrended He II 
intensity data, we do not see a peak in the data at 140 minutes, but 
do see a significant increase at 170 minutes. Thus the He II data 
only weakly support our suggestion that the enhancements seen in 
the Ca II data correspond to chromospheric events that result in the 
oscillations observed with EIS in the Fe lines. 

4. DISCUSSION AND CONCLUSIONS 

As we pointed out in ^ a number of investigations have detected 
Doppler shift oscillations with EIS. Based on the phase differences 
between the Doppler shift and the intensity, we believe that the sig- 
nals we have detected are upwardly propagating magnetoacoustic 
waves. The periods we detect are between 7 and 14 minutes. For 
a set of observations that begins at a particular time, there is con- 
siderable scatter in the measured periods and amplitudes. This is 
probably due to the relatively weak signal we are analyzing. But it 
may also be an indication that a simple sine wave fit is not a good 
representation of the data. It is likely that each line-of-sight passes 
through a complex, time-dependent dynamical system. While a sin- 
gle flux tube may respond to an oscillatory signal by exhibiting a 
damped sine wave in the Doppler shift, a more complex line-of-sight 
may display a superposition of waves. 

Coalignment of the EIS data with both SOT and MDI magne- 
tograms shows that the portion of the EIS slit analyzed in this study 
corresponds to a unipolar flux concentration. SOT Ca II images show 
that the intensity of this feature exhibits 5 -minute oscillations typical 
of chromospheric plasma, but also exhibits some evidence for longer 
period oscillations in the time range detected by EIS. Moreover, the 
Ca II intensity data show that the oscillations observed in EIS are re- 
lated to significant enhancements in the Ca II intensities, suggesting 
that a small chromospheric heating event triggered the observed EIS 
response^ 

Wan g et al] (1200931) also detected propagating slow magnetoa- 
coustic waves in an active region observed with EIS, which they 
associated with the footpoint of a coronal loop. While the oscilla- 
tion periods they measured — 5 minutes — were smaller than those de- 
tected here, many of the overall characteristics we see are the same. 
In each case, the oscillation only persists for a few cycles and the 
phase relationship indicates an upwardly propagating wave. In con- 
trast with their results, however, we do not see a consistent trend for 
the oscillation amplitude to decrease with increasing temperature of 
line formation. Examination of both the Doppler shift data in Fig- 
ure [6] and the results in Table [21 shows that in one case the amplitude 
has a tendency to decrease with increasing temperature of line for- 
mation (oscillation beginning at 113.9 minutes) and in another case 
the amplitude clearly increases with increasing temperature of line 
formation (oscillation beginning at 190 mi nutes). Thu s it does not 
appear th at the results reported by Wang e t al.l (I2009al) are always 
the case. lO'Shea et al.l 0200 2) noted that for oscillations observed 
above a sunspot the amplitude decreased with increasing tempera- 
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ture until the temperature of formation of Mg X, which is formed at 
roughly 1 MK. They then saw an increase in amplitude in emission 
from Fe XVL All the EIS lines we have included in this study have 
temperatures of formation greater than 1 MK. 

Combined EIS and SUMER polar coronal hole observations 
have also shown evide nce propagating slow magnetoacoustic waves 
dBaneriee et alj [200^ . These waves have periods in the 10-30 
minute range. These waves appear to be more like those we observe 
in that t hey are have periods longer than those studied bv lWang et al.l 
([2009ah. 

" Wang et all ([2009ah suggested that the waves they observed were 
the result of leakage of photospheric p-mode oscillations upward into 
the corona. The longer periods we and Banerjee et al. ( 2009) ob- 
serve are probably not related to p-modes. Instead, we speculate 
that the periods of the waves are related to the impulsive heating 
which may be producing them. If an instability is near the point 
where rather than generating a catastrophic release of energy it wan- 
ders back and forth between generating heating and turning back off, 
waves would be created. The heating source could be at a single lo- 
cation, or, for example, locations near each other where instability in 
one place causes a second nearby location to go unstable and begin 
heating plasma. In this view, the periods provide some insight into 
the timescale for the heating to rise and fall and thus may be able to 
place limits on possible heating mechanisms. 

The behavior of slow magnetoacoustic oscillations as a function of 
temperature has been the subject of considerable theoretical work. It 
is generally believed that t he dampin g of the waves i s due to t hermal 
conduction (e.g., De Moo rtel & Hoo d 2003; Klimc huk et al.l l2004). 
Because thermal conduction scales as a high power of the temper- 
ature, conductive damping should be stronger for oscillations de- 
tected in hi gher temperature emission lines (e.g.. Porter et al.lfT994i : 
lOfman & W ang 2002). Earlier EIS observations of the damping of 
standing slow magne toacoustic waves, however, show that this is not 
always the case (Mar iska et al.| [2008). Thus, the temperature behav- 
ior of both the amplitude of the oscillations and the damping, differ 
from some earlier results. We believe that additional observations 
will be required to understand fully the physical picture of what is 
occurring in the low corona when oscillations are observed. Given 
the complex set of structures that may be in the line of sight to any 
given solar location under the EIS slit, we are not entirely surprised 
that different data sets should yield different results, which in some 
cases differ from models. For example, none of the current models 
for oscillations in the outer layers of the solar atmosphere take into 
account the possibility that what appear to be single structures in 
the data might actually be bundles of threads with differing physical 
conditions. 

Our observations along with others (e.g., Wang et al.l I2009al lbl: 
iBaneriee et al.l20 0^ show that low- amplitude upwardly propagating 
slow magnetoacoustic waves are not uncommon in the low corona. 
The periods observed to date range from 5 minutes to 30 minutes. In 
all cases, however, the wave amplitudes are too small to contribute 
significantly to coronal heating. But understanding how the waves 
are generated and behave as a function of line formation tempera- 
ture and the structure of the magnetic field should lead to a more 
complete understanding of the structure of the low corona and its 
connection with the underlying portions of the atmosphere. Instru- 
ments like those on Hinode that can simultaneously observe both the 
chromosphere and the corona, should provide valuable additional in- 
sight into these waves as the new solar cycle rises and more active 



regions become available for study. 
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